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1 Old stories and new measurements 



Enrico Fermi was born hundred yers ago, on September 29 1901. I then take the chance 
of starting my contribution to the COSMO-01 meeting by reminding that the first spec- 
ulations concerning the existence of large scale magnetic fields are due to Fermi who 
published a single author paper in the issue 75 of Physical Review |l[. It was 1949. At 
that time a debate was going on in the community. On one side Alfven ^ and, si- 
multaneously, Richtmyer and Teller |^ were claiming that high energy cosmic rays are in 
equilibrium with stars [@, ^. In a different perspective Fermi gave arguments supporting 
the idea that high energy cosmic rays are in equilibrium with the galaxy. The idea of 
Fermi was that that cosmic rays are a global galactic phenomenon. The suggestion of 
Alven was that cosmic rays are a local solar phenomenon. In order to make his argument 
consistent Fermi needed a large scale (galactic) magnetic field as large as the /i Gauss 
over the size of the galaxy. Four years later, Fermi and Chandrasekar ^ developed the 
first theory of the gravitational instability in the presence of large scale magnetic fields : 
an ancestor of what we call today dynamo theory. 

Fifty years after Fermi's speculations, large scale magnetic fields represent an intrigu- 
ing triple point where cosmic ray physics, cosmology and astrophysics meet for different 
(but related) purposes. Today large scale magnetic fields are measured not only in the 
Milky Way but also in other members of the local group and even in regular Abell clusters. 

From the experimental point of view the best studied field is the one of the Milky Way 
where various observational techniques can be exploited. In particular Zeeman splitting 
estimates offer a good tool in order to study magnetic fields which are locally strong, 
like, for instance, the ones present in compact OH sources. Unfortunately, due to the 
well known limitations of Doppler broadening, Zeeman splitting estimates fail for large 
portions of the interstellar medium where the hyperfine splitting of neutral hydrogen 
would represent an ideal line in order to infer the strength of the magnetic field. Faraday 
rotation measurements, combined with synchrotron emission, represent probably the best 
option in order to measure large scale magnetic fields even outside our galaxy. 

The known limitation of Faraday rotation stems from the need of an independent 
measure of the electron density along the line of sight: today there is no clear experimental 
evidence of why the magnetic field right outside the galaxy should be orders of magnitude 
smaller than the /i Gauss ( n Gauss, as often speculated). 

Since various models of large scale magnetic field generation predict the existence of 
magnetic fields not only in galaxies but also inside clusters, it would be interesting to know 
if magnetic fields are present inside regular Abell clusters. Various results in this direction 
have been reported |0, 0, H, 3. Some studies during the past decade 0, ^ dealt mainly 
with the case of a single cluster (more specifically the Coma cluster). Many radio sources 
inside the cluster were targeted with Faraday rotation measurements (RM). The study 
of many radio sources inside different clusters presents experimental problems due to the 
sensitivity limitations of radio-astronomical facilities. Hence the strategy is currently to 
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study a sample of clusters each with one or two bright radio- sources inside. 

In the past it was shown that regular clusters have a cores with a detectable component 
of RM H, ||. Recent results suggest that fi Gauss magnetic fields are indeed detected inside 
regular clusters |T^. Inside the cluster means in the inter-cluster medium. Therefore, 
these magnetic fields cannot be associated with individual galaxies. 

Regular Abell clusters with strong x-ray emission were studied using a twofold tech- 
nique [0, |ll|. From the ROSAT full sky survey the electron density has been deter- 
mined. Faraday RM (for the same set of 16 Abell clusters) has been estimated through 
observations at the VLA The amusing result (confirming previous claims based only on 
one cluster |P, 0]) is that x-ray bright Abell clusters possess a magnetic field of /x Gauss 
strength. The clusters have been selected in order to show similar morphological features. 
All the 16 clusters monitored with this technique are at low red-shift {z < 0.1) and at 
high galactic latitude (|&| > 20°). 

These recent developments are rather promising and establish a clear connection be- 
tween radio-astronomical techniques and the improvements in the knowledge of x-ray 
sky. There are various satellite missions mapping the x-ray sky at low energies (ASCA, 
CHANDRA, NEWTON [p. There is the hope that a more precise knowledge of the surface 
brightness of regular clusters will help in the experimental determination of large scale 
magnetic fields between galaxies. 

The present paper is organized as follows. In Section II the main aspects of the 
amplification of large scale magnetic fields in the galaxy will be briefly outlined in a 
critical way. In Section III the attention will be concentrated on the general features of 
cosmological magnetogenesis. It will be argued that causal and inflationary mechanisms 
are not alaternative but complementary. In Section IV and V the interplay between the 
evolution of the gauge couplings and the generation of large scale magnetic fields will 
be reviewed. Section IV is devoted to the case of dynamical extra-dimensions leading 
to an effective evolution of the gauge couplings. Section V deals with a model where 
gauge couplings evolve in four dimensions and in a standard framework of cosmological 
evolution. Section VI addresses the possible effects of primordial magnetic fields on the 
polarization of the Cosmic Microwave Background (CMB). In Section VII it will be shown 
that the existence of primordial magnetic fields can be used in order to generate the baryon 
asymmetry of the Universe (BAU). At the same time, the hypermagnetic fields present at 
the electroweak epoch can lead to a gravitational waves background not only in the LISA 
frequency range, but also at higher frequencies. 

*The ROetgen SATellite was flying from June 1991 to February 1999. ROSAT provided a map of the 
x-ray sky in the range 0.1-2.5 keV. 

tThe Very Large Array telescope is a radio-astronomical facility consisting of 27 parabolic antennas 
spread around 20 km in the New Mexico desert. 

*ASCA is operating between 0.4 AND 10 keV and it is flying since February 1993. CHANDRA (NASA 
mission) and NEWTON (ESA mission) have an energy range comparable with the one of ASCA and were 
launched, almost simultaneously, in 1999. 
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2 The galaxy as a gravitationally bound tokamak 



To address the problem of the origin of large scale magnetic fields means to write down 
the equations describing their evolution. This aspect is relatively well understood at least 
for what concerns the late stages of the evolution of the galaxy. A number of excellent 



textbooks and reviews can be consulted [|1^, 0, |1J, 0, |T^. Here only few important 
points will be made clear. 

As far as the evolution of magnetic fields is concerned, the galaxy is a gravitationally 
bound system formed by fluid of charged particles which is globally neutral for scales 
larger than the Debye sphere. In the interstellar medium, where the electron density is 
approximately 3 x 10^^ cm~'^ (as it can be estimated from the dispersion of pulsar signals) 
the Debye sphere has a radius of roughly 10 m. 

Moreover, the galaxy is rotating with a typical rotation period of 3 x 10^ yrs. The 
evolution equations of this system are, physically, the same equations describing the dy- 
namics of electromagnetic fields inside a tokamak. As in the case of a tokamak we have 



two choices. We can study the full kinetic system (the Vlasov-landau equations |TB 



or we can rely on the magnetohydrodynamical (MHD) treatment. Already in flat space 
20| , and, a fortiori, in curved space ||2l]], the kinetic approach is important once we deal 



with electric fields dissipation, charge and current density fluctuations and, in more gen- 
eral terms, with all the high frequency and small length scale phenomena in the plasma 
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Consider a conformally flat Friedmann-Robertson- Walker (FRW) metric written using 
the conformal time coordinate 

ds'^ = a\r])[dT]'^ - dx% (2.1) 

Furthermore, consider an equilibrium homogeneous and isotropic conducting plasma, 
characterized by a distribution function fo{p) common for both positively and negatively 
charged ultrarelativistic particles (for example, electrons and positrons) . Suppose now 
that this plasma is slightly perturbed, so that the distribution functions are 

U{x,p,r]) = fo{p) + 6U{x,p,r]), f_{x,p,r]) = fo{p) + 6f-{x,p,r]), (2.2) 

where + refers to positrons and — to electrons, and p is the conformal momentum. The 
Vlasov equation defining the curved-space evolution of the perturbed distributions can be 
written as [|l^ 

^ + v-^ + e{E + vx B)-^ = (2 3) 

dr] dx dp \ drj /cou' 

^ + ^.^^e{E + vxB)-^=(^) (2 4) 

dr] dx dp \ dr] /coii' 

where the two terms appearing at the right hand side of each equation are the collision 
terms. The electric and magnetic fields are rescaled by the second power of the scale 
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factor. This system of equation represents the curved space extension of the Vlasov- 



Landau approach to plasma fluctuations [0, |T9|. All particle number densities here are 
related to the comoving volume. By subtracting Eqs. ( p. 3D and (|2.4| ) we obtain the 
equations relating the fluctuations of the distributions functions of the charged particles 
present in the plasma to the induced gauge field fluctuations: 

d d ^ dfo 

-Q^fix,p,t) + v- —f{x,p,t) + 2eE--^ = -iy{p)f, 

V ■ E = e j d^pf{x,p,r]), 

V X E + B' = 0, 
V-5 = 0, 

V X B-E' = J £pvf{x,p,ri), (2.5) 
where f{x,p,ri) = 6f+{x,p,ri) — 6f^{x,p,ri) and u^p) is a typical frequency of collisions 
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Now, if Sf±{x,p,ri) 7^ at the beginning of the radiation dominated epoch rjo and 
E{x, ?7o) — B{x, rjo) = initially, the magnetic field at later times can be found from Eqs. 
( p.5| ) [pO|] . Various useful generalizations of the Vlasov-Landau system to curved spaces 
is given in 0, ^, 

For scales sufficiently large compared with the Debye sphere and for frequencies suf- 
ficiently small compared with the plasma frequency the spectrum of plasma excitations 
obtained from the kinetic theory matches the MHD spectrum. 

Since the galaxy is rotating and since the conditions of validity of the MHD approx- 
imation are met, it is possible to use the so-called dynamo instability in order amplify a 
small magnetic inhomogeneity up to the observed value. A necessary condition in order to 
implement this idea is that the flow should not be mirror symmetric, i.e. {v -V x v) ^ 0, 
where v is the bulk velocity of the plasma. The suggestion that the mean magnetic field, 
in a randomly moving medium ( with non mirror-symmetric flow), could grow was first 



proposed by Parker [27]. For a recent (critical) review on galactic dynamos see |T( 
MHD equations can be derived from a microscopic (kinetic) approach and also from a 
macroscopic approach where the displacement current is neglected If the displace- 
ment current is neglected the electric field can be expressed using the Ohm law and the 
magnetic diffusivity equation is obtained 

dB 1 

= ^ x{vxB) + -V^B. (2.6) 
or] a 

The conductivity a appearing in Eq. ( p.6| ) is a global quantity which can be computed 
in a kinetic approach during a given phase of evolution of the background geometry 



P^ . The term containing the bulk velocity field is called dynamo term and it receives 
contribution provided parity is globally broken over the physical size of the plasma. In Eq. 
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( p.6|) the contribution containing the conductivity is usually called magnetic diffusivity 



term. 

The ratio of the two terms on the r.h.s. of Eq. ( p.6|) defines the magnetic Reynolds 

number 

crlV X V X B\ 

ReM ^ Lb a. (2.7) 

\V B\ 

If -R <^ 1 (for a given length scale Lb) the flux lines of the magnetic field will diffuse 
through the plasma. If i? 3> 1 the flux lines of the magnetic field will be frozen into 
the plasma element. From the magnetic diffusivity equation (|2.6| ) it is possible to derive 
the typical structure of the dynamo term by carefully averaging over the velocity field 
according to the procedure outlined in [0, |13[ . By assuming that the motion of the fluid 
is random and that it has zero mean velocity, it is possible to average over the ensemble of 
the possible velocity fields. In more physical terms this averaging procedure of Eq. ( |2.6| ) 
is equivalent to averaging over scales and times exceeding the characteristic correlation 
scale and time tq of the velocity field. This procedure assumes that the correlation scale 
of the magnetic field is much larger than the correlation scale of the velocity field. In this 
approximation the magnetic diffusivity equation can be written as: 

f)R 1 

— = a{VxB) + -V^B, (2.8) 
or] a 

(a = —^{v-'Vxv) is the so-called dynamo term, which vanishes in the absence of 
vorticity; in this equation B is the magnetic field averaged over times larger than tq, 
which is the typical correlation time of the velocity field). The crucial requirement for 
the described averaging procedure is that the turbulent velocity field has to be "globally" 
non-mirror-symmetric. It is interesting to point out []T3[ that the dynamo term in Eq. 
( p.8| ) has a simple electrodynamical meaning, namely, it can be interpreted as a mean 
ohmic current directed along the magnetic field: 

/= -aB. (2.9) 

This equation tells us that an ensemble of screw-like vortices with zero mean helicity is able 
to generate loops in the magnetic flux tubes in a plane orthogonal to the one of the original 
field. This observation will be of some related interest for the physical interpretation of 
the results we are going to present in the following paragraph. We finally notice that if the 
velocity field is parity-invariant (i.e. no vorticity for scales comparable with the correlation 
length of the magnetic field), then the dynamics of the infrared modes is decoupled from 
the velocity field since, over those scales, a = 0. When the (averaged) dynamo term 
dominates in Eq. ( |2.8| ), magnetic fields can be exponentially amplified. When the value 
of the magnetic field reaches the equipartition value (i.e. when the magnetic and kinetic 
energy of the plasma are comparable), the dynamo "saturates". The precise way in which 



the dynamo effect stops was a subject of active research some years ago |12]. What 
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happens is that close to equipartition, Eq. (|2.8| ) should be supplemented with non-linear 
terms whose effect is to stabilize the amplification of the magnetic field to a constant 
value m. 

MHD equations can be studied in two different limits : the ideal (or superconducting) 
approximation where the conductivity is assumed to be very high and the real (or resistive) 
limit where the conductivity takes a finite value. In the ideal limit both the magnetic flux 
and the magnetic helicity are conserved. This means, formally [p^] . 



^ ^ 5 ■ = -- / V X V X 5 ■ dS, (2.10) 



dr] Jt. a Jt. 

where S is an arbitrary closed surface which moves with the plasma. 

If we are in the in the inertial regime (i.e. L > L„ where L„ is the magnetic diffusivity 
length) we can say that the expression appearing at the right hand side is sub-leading and 
the magnetic flux lines evolve glued to the plasma element. Today the magnetic diffusivity 
scale, estimated from MHD considerations, is of the order of the A. U.. This means that 
fields cohernt over distances smaller than the A. U. are dissipated. Conversely, if the 
coherence scale of the magnetic fields is larger than the A. U. the associated magnetic 
flux is conserved. 

It is now worth mentioning the magnetic helicity 



Hm = / d^xA -B, (2.11) 

V 



where A is the vector potential |^ In Eq. ( |2.11| ) the vector potential appears and, therefore 



it might seem that the expression is not gauge invariant. This is not the case. In fact 
A-B is not gauge invariant but, none the less, Hm is gauge-invariant since the integration 
volume is defined in such a way that the magnetic field B is parallel to the surface which 
bounds V and which we will call dV. In n is the unit vector normal to dV then B ■ n in 
dV. 

The magnetic gyrotropy 

B-VxB (2.12) 

it is a gauge invariant measure of the diffusion rate of Hm at finite conductivity. In fact 
2§ 

-^Hm = -- [ d^x B-V X B. (2.13) 



dr] a Jv 

The magnetic gyrotropy is a useful quantity in order to distinguish different mechanisms 
for the magnetic field generation. Some mechanisms are able to produce magnetic fields 
whose fiux lines have a topologically non-trivial structure (i.e. {B ■ V x B) 7^ 0). This 
observation will be used in the following Sections. 

§ Notice that in conformally flat FRW spaces the radiation gauge is conformany invariant. This 
property is not shared by the Lorentz gauge condition p9||. 



6 



The discovery of large scale magnetic fields in the intra-cluster medium implies some 
interesting problems for the mechanisms of generation of large scale magnetic fields. Con- 
sider, first, magnetic fields in galaxies. Usually the picture for the formation of galactic 
magnetic fields is related to the possibility of implementing the dynamo mechanism. By 
comparing the rotation period with the age of the galaxy (for a Universe with Q/^ ~ 0.7, 
h ~ 0.65 and ~ 0.3) the number of rotations performed by the galaxy since its origin 
is approximately 30. During these 30 rotations the dynamo term of Eq. (|2.6| ) dominates 
against the magnetic diffusivity term since parity is globally broken over the physical size 
of the galaxy. As a consequence an instability develops. This instability can be used 
in order to drive the magnetic field from some small initial condition up to its observed 
value. Most of the work in the context of the dynamo theory focuses on reproducing the 
correct features of the magnetic field of our galaxy. For instance one could ask the dynamo 
codes to reproduce the specific ratio between the poloidal and toroidal amplitudes of the 
magnetic field of the Milky Way. 

The achievable amplification produced by the dynamo instability can be at most of 
10^^, i.e. e^°. Thus, if the present value of the galactic magnetic field is 10~^ Gauss, its 
value right after the gravitational collapse of the protogalaxy might have been as small 
as 10~^^ Gauss over a typical scale of 30-100 kpc. 

There is a simple way to relate the value of the magnetic fields right after gravitational 
collapse to the value of the magnetic field right before gravitational collapse. Since the 
gravitational collapse occurs at high conductivity the magnetic flux and the magnetic 
helicity are both conserved. Right before the formation of the galaxy a patch of matter 
of roughly 1 Mpc collapses by gravitational instability. Right before the collapse the 
mean energy density of the patch, stored in matter, is of the order of the critical density 
of the Universe. Right after collapse the mean matter density of the protogalaxy is, 
approximately, six orders of magnitude larger than the critical density. 

Since the physical size of the patch decreases from 1 Mpc to 30 kpc the magnetic 
field increases, because of flux conservation, of a factor (pa/pb)^^^ ~ 10^ where pa and pb 
are, respectively the energy densities right after and right before gravitational collapse. 
Henceforth, the correct initial condition in order to turn on the dynamo instability is 
B ~ 10~^^ Gauss over a scale of 1 Mpc, right before gravitational collapse. 

Since the flux is conserved the ratio between the magnetic energy density, Pb{L, rf) and 
the energy density sitting in radiation, p^{ri) is almost constant and therefore, in terms 
of this quantity (which is only scale dependent but not time dependent), the dynamo 
requirement can be rephrased as 

rB(L) = ^^t^ > 10-=^^ L~lMpc. (2.14) 

If the dynamo is not invoked but the galactic magnetic field directly generated through 
some mechanism the correct value to impose at the onset of gravitational collapse is much 
larger. The possible applications of dynamo mechanism to clusters is still under debate 
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and it seems more problematic [|T0|, [TT|, ^ . The typical scale of the gravitational collapse 
of a cluster is larger (roughly by one order of magnitude) than the scale of gravitational 
collapse of the protogalaxy. Furthermore, the mean mass density within the Abell radius 
( ~ 1.5h-^ Mpc) is roughly 10^ larger than the critical density. Consequently, clusters 
rotate less than galaxies since their origin and the value of rB{L) has to be larger than in 
the case of galaxies. Since the details of the dynamo mechanism applied to clusters are 
not clear, at present, it will be required that rB{L) ^ 10^^'* [for instance rB{L) ~ 10^^^]. 



3 Two classes of mechanisms 

In the context of the ideas illustrated in the previous section, large scale galactic mag- 
netic fields are assumed to be the result of the amplification of a primeval seed. It was 
Harrison [^] who suggested that these seeds might have something to do with cosmology 
in the same way as he suggested that the primordial spectrum of gravitational potential 
fluctuations (i.e. the Harrison- Zeldovich spectrum) might be produced in some primordial 
phase of the evolution of the Universe. Since then, several mechanisms have been invoked 



in order to explain the origin of these seeds [36, 38] and few of them are compatible with 



inflationary evolution. It is not my purpose to review here all the different mechanisms 



which have been proposed and good reviews exist already [p3, 33, EM (see also [35]). A 



very incomplete selection of references is, however, reported [|3^, Here I would 



like to summarize the situation in light of the recent theoretical and experimental devel- 
opments. Consider first of all the case of the galaxy. The mechanisms for magnetic field 
generation can be divided, broadly speaking, into two categories: astrophysical []13|, |1| 



and cosmological. The cosmological mechanisms can be divided, in their turn, into causal 
mechanisms (where the magnetic seeds are produced at a given time inside the horizon) 
and inflationary mechanisms where correlations in the magnetic field are produced out- 
side the horizon. Astrophysical mechanisms have always to explain the initial conditions 
of Eq. ( p.6| ). This is because the MHD are linear in the magnetic fields. It is questionable 
if purely astrophysical considerations can set a natural initial condition for the dynamo 
amplification. 



3.1 Turbulence? 

Both classes of mechanisms have problems. Causal mechanisms usually fail in reproducing 
the correct correlation scale of the field whereas inflationary mechanisms have problems 
in reproducing the correct amplitude required in order to turn on successfully the dynamo 
action. In the context of causal mechanisms there are interesting proposals in order to 
enlarge the correlation scale. These proposals have to do with the possible occurrence of 
turbulence in the early Universe. In order to discuss the turbulent features of a magnetized 
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plasma the kinetic Reynolds number 



Re 



V Lr. 



(3.1) 



should be employed together with the magnetic Reynolds number defined in Eq. 
from the relative balance of the dynamo and diffusivity terms. In Eq. u is the thermal 
diffusivity coefficient, v the bulk velocity of the plasma and Ly the typical correlation scale 
of the velocity field. The ratio of the magnetic Reynolds number to the kinetic Reynolds 
number is the Prandtl number 123 



Pr 



M 



Rcm 
'Re 



L 



B 



(3.2) 



Notice that in Eqs. ( ^.7] ) and ( |3.1| ) the correlations scales of the magnetic field and of 
the fluid motion are, in principle, different. Indeed, as discussed in connection with the 
coarse-grained dynamo equation ( |2.8| ), the correlation scale of the fluid motions should 
be, initially, shorter than the magnetic correlation scale. Even assuming that Lb — 
the Prandtl number is typically much larger than one. Consider, for instance, the case of 
the electroweak epoch |3|, |OHT|, ||, H). At this epoch taking R-^ ~ 3 cm we get 
that PrM ~ z/cr ~ 10^ where the bulk velocity of the plasma is of the order of the bubble 
wall velocity at the epoch of the phase transition. 

This means that the early universe is both kinetically and magnetically turbulent. 
The features of magnetic and kinetic turbulence are different. This aspect reflects in a 

in 



If the 



spectrum of fluctuations is different from the usual Kolmogorov spectrum 
Universe is both magnetically and kinetically turbulent it has been speculated that an 
inverse cascade mechanism can occur JO, 41, |5 
in the context of MHD simulations 



This idea was originally put forward 
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The inverse cascade would imply a growth in 
the correlation scale of the magnetic inhomogeneities and it has been shown to occur 
numerically in the approximation of unitary Prandtl number . Speciflc cascade models 
have been also studied p9| , ^ p| , A particularly important role is played, in this 



context, by the initial spectrum of magnetic flelds (the so-called injection spectrum |41|) 
and by the topological properties of the magnetic flux lines. In particular, if the system 
has non vanishing magnetic helicity and magnetic gyrotropy it was suggested that the 
inverse cascade can occur more efficiently 0. Recently simulations have discussed 
the possibility of inverse cascade in realistic MHD models |^ . More analytic discussions 



based on renormalization group approach applied to turbulent MHD seem to be not totally 



consistent with the occurrence of inverse cascade at large scales [H3 . 



3.2 Magnetic fields from vacuum fluctuations 

Large scale magnetic fiuctuations can be generated during the early history of the Universe 
and can go outside the horizon with a mechanism similar to the one required in order to 
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produce fluctuations in the gravitational potential. In this case the correlation scale of 
the magnetic inhomogeneities can be large. However, the typical amplitudes obtainable 
in this class of models may be too small. The key property allowing the amplification 
of the fluctuations of the scalar and tensor modes of the geometry is the fact that the 
corresponding equations of motion are not invariant under Weyl rescaling of a (conformally 
flat) metric of FRW type. In this sense the evolution equations of relic gravitons and of 
the scalar modes of the geometry are not conformally invariant. The evolution equations 
of gauge fields do not share this property. However, if gauge couplings are dynamical, 
there is a natural way of breaking conformal invariance in the evolution equation of gauge 
fields. Interesting examples in this direction are models containing extra-dimensions and 
scalar-tensor theories of gravity where the gauge coupling is, effectively, a scalar degree 
of freedom evolving in a given geometry. 



4 Dynamical extra-dimensions 



The remarkable similarity of the abundances of light elements in different galaxies leads 
to postulate that the Universe had to be dominated by radiation at the moment when 



the light elements were formed, namely for temperatures of approximately 0.1 MeV [JS, 
|50| . Prior to the moment of nucleosynthesis even indirect informations concerning the 
thermodynamical state of our Universe are lacking even if our knowledge of particle physics 
could give us important hints concerning the dynamics of the electroweak phase transition 



m 



The success of big-bang nucleosynthesis (BBN) sets limits on alternative cosmological 
scenarios. Departures from homogeneity and isotropy |Q of the background geometry 
can be successfully constrained. Bounds on the presence of matter-antimatter domains 
of various sizes can be derived 
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55| , p6| . BBN can also set limits on the dynamical 
evolution of internal dimensions |5^, Internal dimensions are an essential ingredient 
of theories attempting the unification of gravitational and gauge interactions in a higher 
dimensional background like Kaluza-Klein theories |^ and superstring theories 
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Defining, respectively, bsBN and bo as the size of the internal dimensions at the BBN 
time and at the present epoch, the maximal variation allowed to the internal scale factor 
from the BBN time can be expressed as hBBN /bo ~ 1 + e where |e| < 10~^ fS^, |5^. The 
bounds on the variation of the internal dimensions during the matter dominated epoch 
are even stronger. Denoting with an over-dot the derivation with respect to the cosmic 
time coordinate, we have that \b/b\ < lO^^ifo where Ho is the present value of the Hubble 
parameter The fact that the time evolution of internal dimensions is so tightly 

constrained for temperatures lower of 1 MeV does not forbid that they could have been 
dynamical prior to that epoch. Moreover, recent observational evidence [^, Q seem 
to imply that the fine structure constant can be changing even today. 

Suppose that prior to BBN internal dimensions were evolving in time and assume, 
for sake of simplicity, that after BBN the internal dimensions have been frozen to their 
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present (constant) value. Consider a homogeneous and anisotropic manifold whose line 
element can be written as 



ds^ = G^iydx^dx" = o?{ri)[drf — 'jijdx^dx^] — b'^{r])'jabdy"'dy'', 
fi,^ = 0, D — 1 = d + n, i,j = l,...,d, a,b = d + 1, ...,d + n. (4.1) 

[rj is the conformal time coordinate related, as usual to the cosmic time t = J a{ri)dri 
; 7jj(x), 'Jabiy) are the metric tensors of two maximally symmetric Euclidean manifolds 
parameterized, respectively, by the "internal" and the "external" coordinates {x*} and 
{y'^}]. The metric of Eq. ([4.1|) describes the situation in which the d external dimensions 
(evolving with scale factor a{ri)) and the n internal ones (evolving with scale factor b{ri)) 



are dynamically decoupled from each other The results of the present investigation 



however, can be easily generalized to the case of n different scale factors in the internal 
manifold. 

Consider now a pure electromagnetic fluctuation decoupled from the sources, repre- 
senting an electromagnetic wave propagating in the ci-dimensional external space such 
that Afj, = Afj,{x,ri), Aa = 0. In the metric given in Eq. ( |4.1| ) the evolution equation of 
the gauge field fluctuations can be written as 



(^v^G"'^G^"F„^) = 0, (4.2) 



where Fa/3 = V[q,A^] is the gauge field strength and G is the determinant of the D di- 
mensional metric. Notice that if n = the space-time is isotropic and, therefore, the 
Maxwell's equations can be reduced (by trivial rescaling) to the fiat space equations. If 
n 7^ we have that the evolution equation of the electromagnetic fluctuations propagat- 
ing in the external ci-dimensional manifold will receive a contribution from the internal 
dimensions which cannot be rescaled away. 

In the radiation gauge {Aq = and VjA* = 0) the evolution the vector potentials can 
be written as 

b' 

+ nTJ^, - V'A = 0, ^=^- (4.3) 

The vector potentials Ai are already rescaled with respect to the (conformally flat) d + 1 
dimensional metric. In terms of the canonical normal modes of oscillations Ai = b^^'^Ai 
the previous equation can be written in a simpler form, namely 

2 

— * Tl Tl 

A': - V{v)A. - VM. = 0, V{v) = -T^ + -T'. (4.4) 

From this set of equations the induced large scale magnetic flelds can be computed in 
various models for the evolution of the internal manifold . 
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5 Dynamical gauge couplings in four-dimensions 



The evolution of the (Abehan) gauge couphng during an inflationary phase of de Sitter 
type drives the growth of the two-point function of the magnetic inhomogeneities |^8|, |6^ . 



The idea that a theory with local gauge-invariance could lead to a consistent variation of 
the Abelian coupling was explored by Dirac and, subsequently, by Teller and Bekenstein 
||70| , [71] , |72| , [7^ . This physical possibility is also realized in superstring-inspired cosmolog- 
ical scenarios where the gauge coupling is related to the expectation value of the dilaton 



field [q^. In the following we propose a model for the evolution of the gauge coupling in 
a standard cosmological scenario where the inflationary phase is followed by a radiation 
and a matter dominated epoch. 

Suppose that a minimally coupled (massive) scalar field (f) evolves in a conformally 
fiat metric of FRW type (|2.1|) . The field is not the inflaton but it evolves during 
different cosmological epochs parametrized by a different form of a{ri). Typically the 
Universe evolves from an inflationary phase of de Sitter (or quasi-de Sitter) type towards 
a radiation dominated phase which is finally replaced by a matter dominated epoch. 

The evolution equation of in the background given by Eq. ([4.1|) can be written as 



a' 



0" + 27^0' + m^a^cj) = 0, H = -, (5.1) 

a 

where Ti is the Hubble factor in conformal time related to the Hubble parameter in cosmic 
time H = a/a as, Ha = Ti (the dot denotes derivative with respect to cosmic time). 
If evolves during an inflationary phase of de Sitter type the scale factor will be 

«(^) = (--)~\ r,<-Vu (5.2) 

where —rji marks the end of the inflationary phase. If m^a^ ^ Ti (i.e. m <^ H) during 
inflation, according to Eq. (p.!]), relaxes as ~ {—v/Vi)^ V < ~Vi- 
Suppose, as an example, that is coupled to an (Abelian) gauge field 

S ^ J d^xV^(f)^F^^F^'r (5.3) 

The normal modes of the hypermagnetic field fluctuations Bi{x,ri) are 

bi{x, 7]) = (l){r])Bi{x, 77), (5.4) 
and their correlation function during the de Sitter phase can then be written as 

V) = I 7WZ^,P^Ak)Hk, r/)6*(fcr/)e^^^ (5.5) 

where 

P..= /^A..-: 

A;2 



(27r)3 
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The normal modes b{k,ri) will evolve as 

b" + - j]b = 0. (5.7) 

Using now the fact that (f) yf the correlation function, during the de Sitter phase grows 
as 

e?,,(f,r/)~(-r/)-^ (5.8) 
for — i> 0^ (i.e. t oo). Thus, gauge field fluctuations grow during the Sitter stage. 



Furthermore, from Eq. (|5.3| ) the magnetic energy density pB{r,ri) [related to the trace 
of Qij{r,r])] also increases for 77 — >■ 0~. Consequently, since the magnetic energy density 
can be amplified during a de Sitter-like stage of expansion, large scale gauge fluctuations 
pushed outside of the horizon can generate the galactic magnetic field. 

5.1 Evolution of the gauge coupling 

The only gauge couphng free to evolve, in the present discussion, is the one associated 
with the hypercharge field leading, after symmetry breaking, to the time variation of 
the electron charge. In a relativistic plasma the conductivity goes, approximately, as 
T/aem where aem is the fine structure constant . If depends on time also the well 
known magnetohydrodynamical equations (MHD). will have to be generalized, leading, 
ultimately, to different mechanism for the relaxation of the magnetic fields. 

If the evolution of the Abelian coupling is parametrized through the minimally coupled 
scalar field 0, the possible constraints pertaining to the evolution of (p are translated into 
constraints on the evolution of the gauge coupling. Massless scalars cannot exist in the 
Universe: they lead to long range forces whose effect should appear in sub-millimiter tests 
of Newton's law. Consequently, the scalar mass should be, at least, larger than 10"'^ eV 



otherwise it would be already excluded W^. Massive scalars are severely constrained 



from cosmology []76| , [77| , [7^] . When the scalar mass is comparable with the Hubble rate 
(i.e. m ~ H) the field starts oscillating coherently with Planckian amplitude and (p decays 
too late big-bang nucleosynthesis (BBN) can be spoiled. 

Initial conditions for the evolution of (p are given during a de Sitter stage of expansion. 
Thus, the homogeneous evolution of (p can be written as 

(piiv) - h - M-f, v<-vu (5.9) 

where (pi is the asymptotic value of (p which may or may not coincide with the minimum 
of V{(p)] (p2 is also an integration constant. Without fine-tuning (pi and (p2 both coincide 
with Mp. During the infiationary phase m <^ Tii/ai = Hi where. Hi < 10~^ Mp is the 
curvature scale at the end of infiation. 

After rji the Universe enters a phase of radiation dominated evolution (possibly pre- 
ceded by a reheating phase) where the curvature scale decreases. When H^ ~ m the 
scalar field starts oscillating coherently with amplitude (pi. 
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During reheating the scale factor evolves as 0(77) ~ 1]°' so that, in this phase, relaxes 

as 

0rh + (^)'"'". Vl<V< Vr, (5.10) 

where rjj- marks the beginning of the radiation dominated phase occurring at a scale 
ifr > "tTT-- In the case of matter-dominated equation of state during reheating a ~ 2. 

For 1] > rjr, the evolution of the field can be exactly solved (in cosmic time) in terms 
of Bessel functions 



0(t) ~ a-^/^{t)V^[AYi/4{mt) + BJi/^imt)], (5.11) 



where A and B are two integration constants. From Eq. ( p. 11 ), ~ constant + t] ^ for 



H > m, and it oscillates for H < m. When H < m the coherent oscillations of start 
and their energy density decreases as a~^. The curvature scale H^. marks the time at 
which the energy density stored in the coherent oscillations equal the energy density of 
the radiation background, namely 



H^MU-)'-m'^l{-f, (5.12) 



ar 



where r/m correspond to the times at which H ^ m. From Eq. ( ]5.12| ) 

H^^i^^Mp. (5.13) 

where ip = 0i/Mp and ^ = m/Mp. The phase of dominance of coherent oscillation ends 
with the decay of at a scale dictated by the strength of gravitational interactions and 
by the mass m, namely 

H^^^Mp. (5.14) 

In order not to spoil the light elements abundances we have to require that iY^ > i^ns 
implying that m > 10 TeV. 

In order not wash-out the baryon asymmetry produced at the electroweak time by 
overproduction of entropy |7^, H^i, > H^^ may be imposed. Since ifgw ~ \/N^T^^/ Mp 
[where N^a = 106.75 is the effective number of (spin) degrees of freedom at ~ 100 
GeV] we obtain m > 10^ TeV. Notice, incidentally, that the time variation of the gauge 
couplings during the electroweak epoch has not been analyzed and it may be relevant in 
order to produce inhomogeneities at the onset of BBN. 

The inhomogeneous modes of should also be taken into account since we have to 
check that further constraints are not introduced. In order to find how many quanta of 
the field are produced by passing from the inflationary phase to a radiation dominated 
phase let us look at the sudden approximation for the transition of a{ri) |8l|. Consider 
the first order fluctuations of the field ip 

ip{x,T]) = ip{r]) +6ip{x,ri), (5.15) 
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whose evolution equation is, in Fourier space, 



< + 27^V'+ [fc +^V]V' = 0, (5.16) 

where ilj{k,ri) is the Fourier component of 6(f{x,r]). 

In the hmit krji <^ 1 the mean number of quanta created by parametric amphfication 
of vacuum fluctuations is 

1/2 



n{k) ^ |c_(fc)|2 = g|A;r/i|-2^(^— j (5.17) 



where g is a numerical coefficient of the order of 10 ^. the energy density of the created 
(massive) quanta can be estimated from 

dp^ = j^—^mn{k), (5.18) 

where u = k/a is the physical momentum. In the case of a de Sitter phase (A = 3/2) the 
typical energy density of the produced fluctuations is 



p^{r^)^qmHt[j^J (^^j (5.19) 

Also the massive fluctuations may become dominant and we have to make sure that they 
become dominant after 99 already decayed. Deflne as the scale at which the massive 
fluctuations become dominant with respect to the radiation background. The scale if* 
can be determined by requiring that p^iji^) ~ p-yiv*) implying that 

which translates into 

H,-q^e^Mp, (5.21) 

where e = Hi/Mp. In order to make sure that the non-relativistic modes will become 
dominant after already decayed we have to impose that if* < H^p which means that 
m > 10^ TeV for iii ~ 10~^ Mp and which is less restrictive than the other constraints 
previously derived in this paper. 

5.2 Magnetogenesis 

The full action describing the problem of the evolution of the gauge coupling in this 
simplifled scenario is 



S = Jd^x 



(5.22) 
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Using Eq. ([4.1|) the equations of motion become 



8V 



— = -VxE 
or] 

d r 



dr] 

f-B = 0. 



+ J = f{(j))VxB, 
J = a{E + vx B) 



(5.23) 

(5.24) 

(5.25) 
(5.26) 



{B = a^B, E = a^S] J = a^j; a = a^a; B, S, j, Oc are the flat-space quantities whereas 
-B, E^ J , a are the curved-space ones; v is the bulk velocity of the plasma). 

In Eqs. ( ^.23|) -( ^.26|) the effect of the conductivity has been included. The current 
density [present in Eq. ( p.23| ) with a term {dja/ dcfyA""] has been eliminated by the 
using Maxwell's equations. During the inflationary phase, for rj < —rji, the role of the 
conductivity shall be neglected. In this case the evolution equation for the canonical 
normal modes of the magnetic field can be derived from the curl of Eq. ( ^.25|) with the 
useofEq. ( ^T^ : 



b" - 



0, 



(5.27) 



where b = y/jB. For rj > —rji the effect of the conductivity is essential. Therefore, the 
correct equations obeyed by the magnetic field will be the generalization of the MHD 
equations whose derivation will be now outlined. 

MHD equations represent an effective description of the plasma dynamics for large 
length scales (compared to the Debye radius) and short frequencies compared to the 
plasma frequency. MHD can be derived from the kinetic (Vlasov-Landau) equations and 
the MHD spectrum indeed reproduces the plasma spectrum up to the Alvfen frequency. 
MHD can be also derived by neglecting the displacement currents in Eq. ( |5.25|) : 



fV xB = J + f'E. 
By now using the Ohm law together with the Bianchi identity we get to 

(1 + 4)^' = V x{vxB) + -\/^B 
sf s 



(5.28) 



(5.29) 



which is the generalization of MHD equations to the case of evolving gauge coupling. The 
quantity s = a/f is constant. The reason for this statement is the following. The rescaled 
conductivity, 

T 

a = a^a = , (5.30) 
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where ttem ~ /^^- Therefore a/ f = s with these rescahngs, is constant. Taking now the 
Fourier transform of the fields appearing in Eq. ( |5.29| ) the solution, for the Fourier modes, 
will be 

B,{k,ri) = B,{k,7]^)e-^ (5.31) 
Consider now, as an example. 



m = (^J^f- (5-32) 

For 7] < —rji the solution of the evolution equation of the magnetic fluctuations is, from 
Eq. (^) , 



2 



where has been chosen in such a way that h{k^r]) \Jk/2e ^^"^ for 77 — > —00. Using 
Eq. z/ = 5/2. 

For f]i < f] < rjr the Universe is reheating. During this phase the conductivity is not 
yet dominant and the fastest growing solution outside the horizon is given, in the case of 
Eq. (PI), by 



VI 



f 



where we assumed, for concreteness, that a = 2 in Eq. ( p.lO| ). For r] > r]^ Eqs. ( |5.29| ) 
should be used. 

The typical present frequency corresponding to the end of the inflationary phase is 
given, at the present time rjQ, by 

c^i (r^o) ~ 10-^ Tdec C^^'^' (5-35) 

Notice that 10~^Tdec = 100 GHz where Tdec is the decoupling temperature. In Eq. (|5.35|) , 
( = Hj./Mp where H,. is the typical curvature scale associated with t],. [see e.g. Eqs. 
( p.ll| )-( |5.12| )]. In the case of instantaneous reheating Hi ~ H^. The typical (present) 
frequency corresponding to the onset of the radiation dominated phase is given by 

'"^^^^^ - (5.36) 



uJiiVo) Ve 

For t] > T]j. the conductivity dominates the evolution and using Eq. (|5.29|) we can 
estimate the trace of the two-point function ( ^l5|) 

pB{r, V) = J PB{k V>—j^, — (5-37) 

with ^ 

PB{k,7])='^\b{k,r])f. (5.38) 
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Thus, in terms of 



and using Eqs.(p:29D-(ra) 



rB{u;,r]o)=C{u;,,u;r)r{u) (5.40) 



where 



c^(-i,c..) = c4r^(5/2)(^)''"'' 



n.) = e-^^[i^)i^)H^f^]^, (5.41) 

UJl LOya UJc 



with cjo- ~ ys/rjo and Tq ~ 10 GeV. In Eq. ( p.41| ) lu^, uj^n and are, respectively, the 
present values of H^, and He- Using the notation of Eq. (|5.35|) we have that 



All the frequencies are evaluated at the time rjQ. Using now the previous equations, 

rsicoG) ~ et"'- (5.43) 

where ujq ~ 10~^^Hz is the present frequency corresponding to a Mpc scale. This result 
should be confronted with typical values of tb required in order to explain galactic (and 
possibly inter-galactic) magnetic fields. Before doing so let us discuss the main spectral 
features of the obtained results. The spectrum is flat for frequencies u < Ua- where u^j is 
the present value of the magnetic diffusivity frequency. The frequency Ua roughly of the 
order of 10"'^ Hz corresponding to a typical scale of the order of the astronomical unit. 
Thanks to the flatness of the spectrum, the constraints obtained on the parameters at the 
galactic scale will be preserved at even larger scales. The upper limit in the amplitude 
will be dictated by e = 10~^ (as required by CMB observations) and ( — e (i.e. the case 
of instantaneous reheating). These features are summarized in Fig. |l| 

In Fig. 1^ the shaded area illustrates the region where magnetogenesis is possible. 

To be consistent with inflationary production of scalar and tensor fluctuations of the 
geometry e < 10~^ should be imposed. Thus, in Fig. ^ the parameters should all lie 
below the (horizontal) thick line. Moreover, since H,. < Hi, ( < e. Recall that ( > 10~^^ 
in order not to affect the nucleosynthesis epoch [see the dot-dashed line in Fig. ||. This 
requirement comes about since > m and ^ > 10^^^. In these considerations have 
been extended to more general evolutions of the gauge coupling. 
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log ( CD /Hz) 



Figure 1: Different spectra are reported for various choices of the parameters. The value 
of ^ has been fixed to 10"^^ in order to have m > lO^TeV, as required by electroweak 
baryogenesis considerations. 



6 Faraday rotation of CMB 

If large scale magnetic fields are truly primordial they should be present prior to the 
decoupling epoch. Can these fields be "detected" in some way? The effects of large 
scale magnetic fields possibly present to the decoupling epoch are manifold. Large scale 
magnetic fields possibly present at the BBN epoch can have an impact on the light nuclei 
formation. By reversing the argument, the success of BBN can be used in order to bound 
the magnetic energy density possibly present at the time of formation of light nuclei. 
These bounds are qualitatively different from the ones previously quoted and coming, 



alternatively, from homogeneity [p2] and isotropy [53| of the background geometry at the 



BBN time. As elaborated in slightly different frameworks through the years 



85| , |86| , magnetic fields possibly present at the BBN epoch could have a twofold effect. On 
one hand they could enhance the rate of reactions (with an effect proportional to aps) 
and, on the other hand they could artificially increase the expansion rate (with an effect 
proportional to ps)- It turns out that the latter effect is probably the most relevant p5[] . 
In order to prevent the Universe from expanding too fast at the BBN epoch ps < 0.27pu 
where pi^ is the energy density contributed by the standard three neutrinos for T < 1 
MeV. 

If magnetic fields are present at BBN there are no reasons why they should be absent 
at the decoupling epoch. Moreover, if quantum mechanical fluctuations of gauge fields 
are amplified by breaking Weyl invariance, then magnetic fields will be produced over 
different scales including scales larger than the Mpc. 

Large scale magnetic fields present at the decoupling epoch can have various con- 
sequences. For instance they can induce fluctuations in the CMB p5| , they can 
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-14 



-12 



Figure 2: Magnetogenesis requirements are illustrated. The thin line b corresponds to 
'"b('^g) ~ 10"'^^. The dashed line a corresponds r^{ujQ) ~ 10~^°. The (horizontal) thick 
line correspond to e = 10~^. The (diagonal) thick line corresponds to e ~ The thick line 
c corresponds to the critical density bound which implies r{uj) < 1. Since the spectrum is 
flat (or mildly increasing with the frequency), this bound is also satisfied at large scales. 
Magnetogenesis is possible when the parameters lie in the shaded area 



distort the Planckian spectrum of CMB they can distort the acoustic peaks of CMB 
anisotropics |^| and they can also depolarize CMB 
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The polarization of the CMB represents a very interesting observable which has been 



extensively investigated in the past both from the theoretical |]91]] and experimental points 
Forthcoming satellite missions like PLANCK |^ seem to be able to achieve 



of view 92 



a level of sensitivity which will enrich decisively our experimental knowledge of the CMB 
polarization with new direct measurements. 

If the background geometry of the universe is homogeneous but not isotropic the 
CMB is naturally polarized ||91|. This phenomenon occurs, for example, in Bianchi-type I 



models @]. On the other hand if the background geometry is homogeneous and isotropic 
(like in the Friedmann-Robertson- Walker [FRW] case) it seems very reasonable that the 
CMB acquires a small degree of linear polarization provided the radiation field has a 
non- vanishing quadrupole component at the moment of last scattering P3 . 

Before decoupling photons, baryons and electrons form a unique fluid which possesses 
only monopole and dipole moments, but not quadrupole. Needless to say, in a homo- 
geneous and isotropic model of FRW type a possible source of linear polarization for 
the CMB becomes efficient only at the decoupling and therefore a small degree of linear 
polarization seems a firmly established theoretical option which will be (hopefully) sub- 
jected to direct tests in the near future. The discovery of a linearly polarized CMB could 
also have a remarkable impact upon other (and related) areas of cosmology. Indeed the 
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linear polarization of the CMB is a very promising laboratory in order to directly probe 
the speculated existence of a large scale magnetic field (coherent over the horizon size at 
the decoupling) which might actually rotate (through the Faraday effect fl^, |T3|, |14|) the 
polarization plane of the CMB. 

Consider, for instance, a linearly polarized electromagnetic wave of physical frequency 
uj traveling along the x direction in a cold plasma of ions and electrons together with a 
magnetic field {B) oriented along an arbitrary direction ( which might coincide with x in 
the simplest case). If we let the polarization vector at the origin (x = ?/ = z = 0, t = 0) 
be directed along the y axis, after the wave has traveled a length Ax, the corresponding 
angular shift (Aa) in the polarization plane will be : 

(conventions: ujb = eB/m is the Larmor frequency; ujpi = Aimee^ / m is the plasma 
frequency is the electron density and fe is the ionization fraction ; we use everywhere 
natural units h = c = ks = It is worth mentioning that the previous estimate of 
the Faraday rotation angle Aa holds provided ^ cu^ and u ^ Upi. From Eq. ( |6.1| ) 
by stochastically averaging over all the possible orientations of B and by assuming that 
the last scattering surface is infinitely thin (i.e. that Axfeiie — cr^^ where ax is the 
Thompson cross section) we get an expression connecting the RMS of the rotation angle 
to the magnitude of i? at t ~ tdec 



{(Aa)Y'^^-(^'\^^][—), 5c = 10-=^ Gauss, um ^SxlO'' Hz (6.2) 




(in the previous equation we implicitly assumed that the frequency of the incident electro- 
magnetic radiation is centered around the maximum of the CMB). We can easily argue 
from Eq. ( |6.2|) that if B{tdec)^Bc the expected rotation in the polarization plane of the 
CMB is non negligible. Even if we are not interested, at this level, in a precise estimate 
of Aa, we point out that more refined determinations of the expected Faraday rotation 
signal (for an incident frequency um ~ 30 GHz) were recently carried out ||96|, ^] leading 
to a result fairly consistent with (|6TT|). 

Then, the statement is the following, //the CMB is linearly polarized and if a large 
scale magnetic field is present at the decoupling epoch, then the polarization plane of the 
CMB can be rotated [Q. The predictions of different models predicting the generation 
of large scale magnetic fields can then be confronted with the requirements coming from 



a possible detection of depolarization of the CMB |90]. 



7 Hypermagnetic fields, EWPT and relic gravitons 

If magnetic fields are generated over all physical scales compatible with the plasma dy- 
namics they may have been present even before the BBN epoch, namely at the electroweak 
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epoch. Some of these fields already decayed by today since they are washed out by simul- 
taneous effects of finite magnetic and thermal diffusivities. However, at the electroweak 
time these fields did not dissipate yet. When we talk about large scale magnetic fields 
at the electroweak scale we use, as unit, the electroweak time at the epoch of the phase 
transition. 

In this section the we will illustrate the idea that a classical hypermagnetic back- 
ground can provide an explanation of the possible formation of the BAU giving rise, 
simultaneously, to stochastic gravitational waves background. 

At small temperatures and small densities of different fermionic charges the SUl{'2) ® 
is broken down to the f/em(l) and the long range fields which can survive in the 
plasma are the ordinary magnetic fields. However, for sufficiently high temperatures the 
SUl{'2.) ® f/y(l) is restored and non-screened vector modes correspond to hypermagnetic 
fields. At the electroweak epoch the typical size of the horizon is of the order of 3 cm . 
The typical diffusion scale is of the order of 10~® cm. Therefore, over roughly eight orders 
of magnitude hypermagnetic fields can be present in the plasma without being dissipated 
p3| . The evolution of hypermagnetic fields can be obtained from the anomalous magne- 
tohydrodynamical (AMHD) equations. The AMHD equations generalize the treatment 
of plasma effects involving hypermagnetic fields to the case of finite fermionic densityp4]. 

Depending on their topology, hypermagnetic fields can have various consequences [43, 
If the hypermagnetic fiux lines have a trivial topology they can have an impact 
on the phase diagram of the electroweak phase transition [|104| , |105|| . If the topology of 
hypermagnetic fields is non trivial, hypermagnetic knots can be formed |^ and, under 
specific conditions, the BAU can be generated p^ . 

A classical hypermagnetic background in the symmetric phase of the EW theory can 
produce interesting amounts of gravitational radiation in a frequency range between 10~"^ 
Hz and the kHz. The lower tail falls into the LISA window while the higher tail falls in 
the VIRGO/LIGO window. For the hypermagnetic background required in order to seed 
the BAU the amplitude of the obtained GW can be even six orders of magnitude larger 
than the infiationary predictions. In this context, the mechanism of baryon asymmetry 
generation is connected with GW production pSl, p9 . 



7.1 Hypermagnetic knots 

It is possible to construct hypermagnetic knot configurations with finite energy and he- 
licity which are localized in space and within typical distance scale L^. Let us consider in 
fact the following configuration in spherical coordinates |91 



2Bn cos 9 



TtLs [7^2 + 1] 

2Bq sin^ 



2 ' 
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(7.1) 



where TZ = r / Lg is the rescaled radius and Bq is some dimensionless amphtude and n is 
just an integer number whose physical interpretation will become clear in a moment. The 
hypermagnetic field can be easily computed from the previous expression and it is 



4i?o cos 9 



2 ' 



VT LI [7^2 + 1] 
8fio 7^sin0 



Tnsin6', 



(7.2) 



The poloidal and toroidal components of 7i can be usefully expressed as Tip = Hrer+Hgeg 
and Ht = 'H(j,es. The Chern- Simons number is finite and it is given by 



9 



12 



CS 



327r2 Jv ^ ■ ^^"^ ^ 327r2 io vr^ [7^2 + if 327r2 



(7.3) 



We can also compute the total helicity of the configuration namely 



nl ^ n-^, ,3 256 Bl n 
V ttL^ jo (l + Tc"')^ 



bB^n 



(7.4) 



We can compute also the total energy of the field 

1 



E 



2Jv ' ' 2 L 



(7.5) 



and we discover that it is proportional to n^. This means that one way of increasing the 
total energy of the field is to increase the number of knots and twists in the fiux lines. We 
can also have some real space pictures of the core of the knot (i.e. TZ = r/Lg < 1). This 
type of configurations can be also obtained by projecting a non-Abelian SU(2) (vacuum) 
gauge field on a fixed electromagnetic direction ||10(]|| [] These configurations have been also 
studied in [ p.02| , lUB ]. In particular, in [ I03|| , the relaxation of HK has been investigated 
with a technique different from the one employed in pSl, 091 but with similar results. 



Topologically non-trivial configurations of the hypermagnetic fiux lines lead to the 
formation of hypermagnetic knots (HK) whose decay might seed the Baryon Asymmetry 

^ In order to interpret these solutions it is very iteresting to make use of the Clebsh decomposition. 
The imphcations of this decomposition (beyond the hydrodynamical context, where it was originally dis- 
covered) have been recently discussed (see [lOl] and references therein). I thank R. Jackiw for interesting 
discussions about this point. 
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of the Universe (BAU). HK can be dynamically generated provided a topologically trivial 
(i.e. stochastic) distribution of flux lines is already present in the symmetric phase of 
the electroweak (EW) theory ||9^, In spite of the mechanism generating the HK, 
their typical size must exceed the diffusivity length scale. In the minimal standard model 
(MSM) (but not necessarily in its supersymmetric extension) HK are washed out. A 
classical hypermagnetic background in the symmetric phase of the EW theory can produce 
interesting amounts of gravitational radiation. 

The importance of the topological properties of long range (Abelian) hj^ercharge 



magnetic fields has been stressed in the past [p.06| , p.071 , p.08| , p.09|| . In [ p.lCI|| it was argued 
that if the spectrum of hypermagnetic fields is dominated by parity non-invariant Chern- 
Simons (CS) condensates, the BAU could be the result of their decay. Most of the 
mechanisms often invoked for the origin of large scale magnetic fields in the early Universe 
seem to imply the production of topologically trivial (i.e. stochastic) configurations of 



magnetic fields |36, 37, 38 



7.2 Hypermagnetic knots and BAU 

Suppose that the EW plasma is filled, for T > with topologically trivial hypermagnetic 
fields Tiy, which can be physically pictured as a collection of fiux tubes (closed because of 
the transversality of the field lines) evolving independently without breaking or intersect- 
ing with each other. If the field distribution is topologically trivial (i.e. (?iyVx?-^y) = 0) 
parity is a good symmetry of the plasma and the field can be completely homogeneous. 
We name hypermagnetic knots those CS condensates carrying a non vanishing (averaged) 
hypermagnetic helicity (i.e. (TYy ■ V x Tiy) 7^ 0). If (?iy ■ V x Tiy) 7^ parity is broken 
for scales comparable with the size of the HK, the fiux lines are knotted and the field Tiy 
cannot be completely homogeneous. 

In order to seed the BAU a network of HK should be present at high temperatures 



44| , |110|| . In fact for temperatures larger than Tc the fermionic number is stored both 
in HK and in real fermions. For T < Tc, the HK should release real fermions since the 
ordinary magnetic fields (present after EW symmetry breaking) do not carry fermionic 
number. If the EWPT is strongly first order the decay of the HK can offer some seeds for 
the BAU generation ||110|| . This last condition can be met in the minimal supersymmetric 
standard model (MSSM) 1TTT|, |m|, |rg, |Tg 



Under these hypotheses the integration of the f/(l)y anomaly equation ||11(]|| gives the 
CS number density carried by the HK which is in turn related to the density of baryonic 
number Ub for the case of fermionic generations. 

ns, , a' nf{nY-^xnY)MoT , g'^ ^. 

{g' is the f/(l)y coupling and s = (2/45)7r^A''e//T^ is the entropy density; Nejf, at Tc, 
is 106.75 in the MSM; Mq = Mp/lM^/NJf ~ 7.1 x lO^^GeV). In Eq. (^ T is 
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the perturbative rate of the right electron chirality flip processes (i.e. scattering of right 
electrons with the Higgs and gauge bosons and with the top quarks because of their large 
Yukawa coupling) which are the slowest reactions in the plasma and 

is the rate of right electron dilution induced by the presence of a hypermagnetic field. In 
the MSM we have that F < Tj^ ||115|| whereas in the MSSM F can naturally be larger than 



F-^. Unfortunately, in the MSM a hypermagnetic field can modify the phase diagram of 
the phase transition but cannot make the phase transition strongly first order for large 
masses of the Higgs boson ||104[. Therefore, we will concentrate on the case F > Ft^ 



and we will show that in the opposite limit the BAU will be anyway small even if some 
(presently unknown) mechanism would make the EWPT strongly first order in the MSM. 

HK can be dynamically generated [0, ^ . Gauge-invariance and transversality of the 
magnetic fields suggest that perhaps the only way of producing {Tiy ■ V x Tiy) 7^ is to 
postulate, a time-dependent interaction between the two (physical) polarizations of the 
hypercharge field Ya. Having defined the Abelian field strength Yap = VfaYg] and its dual 
Yo,f3 such an interaction can be described, in curved space, by the Lagrangian density 



--Y aY"f^ + C—Y 
4 4M "'^ 



(7i 



where Qfj^^y is the metric tensor and g its determinant, c is the coupling constant and M 
is a typical scale. This interaction is plausible if the U{1)y anomaly is coupled, (in the 
symmetric phase of the EW theory ) to dynamical pseudoscalar particles ip (like the axial 
Higgs of the MSSM). Thanks to the presence of pseudoscalar particles, the two polar- 
izations of Tiy evolve in a slightly different way producing, ultimately, inhomogeneous 
HK. 

Suppose that an inflationary phase with a(r) ~ is continuously matched, at the 
transition time ri, to a radiation dominated phase where a(r) ~ r. Consider then a 
massive pseudoscalar field ijj which oscillates during the last stages of the inflationary 
evolution with typical amplitude ipo ~ M. As a result of the inflationary evolution 
iVipl -C ip'. Consequently, the phase of ip can get frozen. Provided the pseudoscalar 
mass m is larger than the inflationary curvature scale Hi ~ const., the ip oscillations are 
converted, at the end of the quasi-de Sitter stage, in a net helicity arising as a result of 
the different evolution of the two (circularly polarized) vector potentials 

Y± + + ^lY± = 0, Hy = VxY (7.9) 
col = k"^ ^ k—aip (7.10) 

(where we denoted with Hy = a'^Hy the curved space fields and with a = the rescaled 
hyperconductivity; the prime denotes derivation with respect to conformal time r whereas 
the over-dot denotes differentiation with respect to cosmic time t). 
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Since ^ uo^ the helicity gets amplified according to Eq. (|7.9|) . There are two 
important points to stress in this context. First of all the plasma effects as well as the 
finite density effects are important. This means, in practical terms, that the dissipation 
scales of the problem should be borne in mind. This has not always been done. The second 
point is related to the first. By treating, consistently, the plasma and finite density effects 
in the context of AMHD |9^, one realizes that the pseudoscalar coupling of if) 
together with the coupling of the coupling to the chemical potential are not sufficient in 
order to seed the BAU unless some hypermagnetic background is already present. In other 
words the scenario which leads to the generation of the BAU is the following. Correlation 
in the hypermagnetic fields are generated outside the horizon durin inflation by direct 
breaking of conformal invariance. An example in this direction has been given in the 
previous Section. Then hypermagnetic fields re-entering at the electroweak epoch will 
participate in the dynamics and, in particular, they will feel the effect of the anomalous 
coupling either to ip or to the chemical potential. The effect of the anomalous coupling 
will not be to amplify the hypermagnetic background. The anomalous coupling will 
make the topology of the hypermagnetic flux lines non-trivial. So the statement is that 
conformal invariance is broken and, hypermagnetic have anomalous couplings, then 
a BAU >10~^'^ can be achieved without spoiling the standard cosmological evolution 
||98| , p9l . It is worth mentioning that this type of scenario may be motivated by the 
low energy string effective action where, by supersymmetry, Kalb-Ramond axions and 
dilatons are coupled, respectively, to the anomaly and to the gauge kinetic term ||116| . 



It is interesting to notice that, in this scenario, the value of the BAU is determined by 
various particle physics parameters but also by the ratio of the hypermagnetic energy 
density over the energy density sitting in radiation during the electroweak epoch, namely. 



using the language of the previous section ^ 

— oc r. (' -ll) 

s 

In order to get a sizable BAU r should be at least 10~^ if the anomalous coupling operates 
during a radiation phase. The value of r coupld be smaller in models where the anomalous 
coupling is relevant during a low scale inflationary phase ||9^ . 



7.3 Hypermagnetic fields and relic gravitons 

If a hypermagnetic background is present for T > T^,, then, as discussed in | |117| | in the 
context of ordinary MHD, the energy momentum tensor will acquire a small anisotropic 
component which will source the evolution equation of the tensor fluctuations h^^, of the 
metric Qf^y. 

hfr + 27^/l^ - V^hij = -IQnGT^p. (7.12) 
where t-j is the tensor component of the energy-momentum tensor of the hyper- 



magnetic fields. Suppose now, as assumed in [104 that has constant amplitude and 
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that it is also homogeneous. Then as argued in ||118|| we can easily deduce the critical 



fraction of energy density present today in relic gravitons of EW origin 

n^{to) = ^^z-y, p,{T,) N^sT^ (7.13) 

Pc 

{zeq = 6000 is the redshift from the time of matter-radiation, equality to the present time 
to)- Because of the structure of the AMHD equations, stable hypermagnetic fields will 
be present not only for cUew ~ Kw/d but for all the range tUew < w < uj„ where Ua- is the 
diffusivity frequancy. Let us assume, for instance, that Tc ~ 100 GeV and N^ff = 106.75. 
Then, the (present) values of a;ew is 



...(to).2.01xlO-^^J^^J Hz. (7.14) 

Thus, co'o-(to) ~ lO^cUew Suppose now that ~ 100 GeV; than we will have that to'ew(to) ~ 
10~^ Hz. Suppose now, as assumed in ||104|| , that 

\n\/T^>0.3. (7.15) 

This requirement imposes r ~ 0.1-0.001 and, consequently, 

hlQcy^ ^ 10"^ - 10-^ (7.16) 

Notice that this signal would occurr in a (present) frequency range between 10~^ and 
10^ Hz. This signal satisfies the presently available phenomenological bounds on the 
graviton backgrounds of primordial origin. The pulsar timing bound ( which applies for 
present frequencies up ~ 10~* Hz and implies /igl^cw < 10~^) is automatically satisfied 
since our hypermagnetic background is defined for 10~^Hz < a; < lO'^Hz. The large 
scale bounds would imply /iq^qw < 7 x 10~^^ but a at much lower frequency (i.e. 10~^^ 
Hz). The signal discussed here is completely absent for frequencies u < uj^^. Notice 
that this signal is clearly distinguishable from other stochastic backgrounds occurring at 
much higher frequencies (GHz region) like the ones predicted by quintessential inflation 
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119|| . It is equally distinguishable from signals due to pre-big-bang cosmology (mainly 



in the window of ground based interferometers [1^]). The frequency of operation of the 



interferometric devices (VIRGO/LIGO) is located between few Hz and 10 kHz ||120|| . The 
frequency of operation of LISA is well below the Hz. Hence, the background described 
here lies in the LISA window |p8| , |99| . In this model the signal can be located both in 
the LISA window and in the VIRGO/LIGO window due to the hierarchy between the 
hypermagnetic diffusivity scale and the horizon scale at the phase transition. 



8 Concluding remarks 

Fifty years after the first seminal speculations on the existence of large scale magnetic 
fields there are still various dark corners. It is crucial to improve the present evidence 



27 



on large scale magnetic fields inside clusters. The existence of // Gauss magnetic fields 
inside clusters would be interesting both for magnetogenesis models.The ongoing progress 
in x-ray astronomy (driven by satellite missions) seems to be rather promising in order to 
improve our knowledge of the density of charged particles inside the clusters. Even more 
important would be to obtain observational evidence of magnetic fields right outside the 
Milky Way. At the moment there is no conclusive evidence of the fact that these fields 
should be as small as the n Gauss. They can be indeed larger. Another reasonable 
question often asked discussing large scale magnetic fields concerns the possibility of 
their primordial origin. In this direction an interesting observable would certainly be the 
polarization of the CMB. Moreover, large scale magnetic fields present prior to BBN could 
give rise to an interesting stochastic GW background in the LISA window. 

The detection of large scale magnetic fields inside clusters puts pressure on many 
models of magnetogenesis. Few models predict magnetic fields of reasonable magnitude 
in this respect. A more complete theoretical understanding of the dynamo mechanism in 
clusters is certainly needed. 

As far as the problem of primordial magnetogenesis is concerned. Few remarks are in 
order. Causal mechanisms may work at least for the galactic magnetic fields (where the 
understanding of the dynamo action is more convincing) provided some inverse cascade 
is realized. As far as inflationary mechanisms are concerned a crucial problem is to find 
a consistent way in order to break conformal invariance. If gauge couplings are dynami- 
cal, conformal invariance may be plausibly broken in the early stages of the inflationary 
evolution. 



ACKNOWLEDGEMENTS 

The author wishes to thank M. E. Shaposhnikov for inspiring collaboration. Interesting 
discussions with M. Laine, H. Kurki-Suonio, E. Sihvola, K. Jedamzik are also acknowl- 
edged. Finally the author wishes to thank K. Enqvist, K. Kajantie, K. Kainulainen and 
M. Ross for the kind invitation to the COSMO-01 meeting. 



28 



References 

[1] E. Fermi, Phys. Rev. 75, 1169 (1949). 

[2] H. Alfven, Arkiv. Mat. F. Astr., o. Fys. 29 B, 2 (1943). 

[3] H. Alfven, Phys. Rev. 75, 1732 (1949). 

[4] R. D. Richtmyer and E. Teller Phys. Rev. 75, 1729 (1949). 

[5] E. Fermi and Chandrasekar Astrophys. J. 118, 113 (1953). 

[6] K.-T. Kim, P.P. Kronberg, P. D. Dewdney and T. L. Landecker, Astrophys. J. 355, 
29 (1990). 

[7] L. Ferretti, D. Dallacasa, G. Giovannini, and A. Tagliani, Astronom. and Astrophys. 
302, 680 (1995). 

[8] K.-T. Kim, P. C. Tribble, and P.P. Kronberg, Astrophys. J. 379, 80 (1991). 

[9] O. Goldshmidt and Y. Raphaeh, Astrophys. J. 411, 518 (1993). 
[10] T.E. Clarke, P.P. Kronberg and H. Bohringer |astro-ph /001128l| . 
[11] H. Bohringer, Rev. Mod. Astron. 8, 295 (1995). 

[12] Ya. B. Zeldovich, A. A. Ruzmaikin, and D.D. Sokoloff, Magnetic Fields in Astro- 
physics (Gordon and Breach Science, New York, 1983). 

[13] E. N. Parker, Cosmical Magnetic Fields (Clarendon, Oxford, 1979). 

[14] P. P. Kronberg, Rep. Prog. Phys. 57, 325 (1994). 

[15] R. Beck, A. Brandenburg, D. Moss, A. Skhurov, and D. Sokoloff Annu. Rev. Astron. 
Astrophys. 34 , 155 (1996). 

[16] R. M. Kulsrud, Annu. Rev. Astron. Astrophys. 37, 37 (1999). 

[17] M. Giovannini and M. Shaposhnikov, Phys. Rev. D 62, 103512 (2000). 

[18] A. Vlasov, Zh. Eksp. Teor. Fiz. 8, 291 (1938); J. Phys. 9, 25 (1945). 

[19] L. D. Landau, J. Phys. U.S.S.R. 10, 25 (1945). 

[20] E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics^ (Pergamon Press, Oxford, 
England, 1980). 



29 



[21] J. Bernstein, Kinetic Theory in the Expanding Universe, (Cambridge University 
Press, Cambridge, England, 1988). 

[22] N. A. Krall and A. W. Trivelpiece, Principles of Plasma Physics, (San Francisco 
Press, San Francisco 1986). 

[23] D. Biskamp, Non-linear Magnetohydrodynamics (Cambridge University Press, Cam- 
bridge, 1994). 

[24] C. P. Dettmann, N. E. Frankel, and V. Kowalenko, Phys. Rev. D 48, 5655 (1993). 

[25] R. M. Gailis, C. P. Dettmann, N. E. Frankel, and V. Kowalenko, Phys. Rev. D 50, 
3847 (1994). 

[26] R. M. Gailis, N. E. Frankel, and C. P. Dettmann, Phys. Rev. D 52, 6901 (1995). 
[27] E. N. Parker, Astrophys. J. 122, 293 (1955). 
[28] M. Giovannini, Phys. Rev. D 58, 12407 (1998). 
[29] L. H. Ford, Phys.Rev. D 31, 704 (1985). 

[30] S. A. Colgate, H. Li, and V. Pariev, |astro-ph/001"248^ ; S. A. Colgate and H. Li, 
|astro-ph/0001418i 

[31] E. R. Harrison, Phys. Rev. Lett. 30, 188 (1973). 

[32] K. Enqvist, Int. J.Mod.Phys.D 7, 331 (1998). 

[33] D. Grasso and H. Rubinstein, Phys. Rept. 348, 163 (2001). 

[34] A. D. Dolgov, |hep-ph/0110293| . 

[35] T. Prokopec |astro-ph/0106247 . 

[36] M. S. Turner and L. M. Widrow, Phys. Rev. D 37, 2734 (1988);T. Vachaspati and 
A. Vilenkin, Phys.Rev.Lett. 67, 1057 (1991); B. Ratra, Astrophys. J. Lett. 391, LI 
(1992); S. Carroll, G. Field and R. Jackiw, Phys. Rev. D 41, 1231 ; W. D. Garretson, 
G. Field and S. Carroll, Phys. Rev. D 46, 5346 (1992); O. Bertolami and D. Mota, 
Phys.Lett.B 455, 96 (1999); A. C. davis and K. Dimopoulos, Phys. Rev. D 55, 7398 
(1997). 

[37] A. Dolgov and J. Silk, Phys. Rev. D 47, 3144 (1993); A. Dolgov, Phys.Rev.D 48, 
2499 (1993). 



30 



[38] K. Enqvist and P. Olesen, Phys. Lett. B 319, 178 (1993); ibtd. 329, 195 (1994); T. 
W. Kibble and A. Vilenkin, Phys. Rev. D 52, 679 (1995); G. Baym, D. Bodeker and 
L. McLerran, ibid. 53, 662 (1996); J. M. Quashnock, A. Loeb and D. N. Spergel, 
Astrophys. J. 344, L49 (1989); B. Cheng and A. V. Ohnto, Phys. Rev. D 50, 2421 
(1994). 

[39] A. Brandenburg, K. Enqvist, and P. Olesen, Phys. Rev. D 54, 1291 (1996). 

[40] A. Brandenburg, K. Enqvist, and P. Olesen, Phys. Lett.B 392, 395 (1997) 

[41] P. Olesen, Phys. Lett. B 398, 321 (1997). 

[42] K. Dimopoulos and A. C. Davis, Phys.Lett. B 390, 87 (1997). 

[43] M. Giovannini and M. Shaposhnikov, Phys. Rev. D 57, 2186 (1998). 

[44] M. Giovannini and M. Shaposhnikov, Phys. Rev. Lett. 80, 22 (1998). 

[45] D. T. Son, Phys. Rev. D 59, 063008 (1999). 

[46] A. Pouquet, U. Frish, and J. Leorat, J. Fluid Mech. 77, 321 (1976); J. Leorat, A. 
Pouquet, and U. Frish ibid. 104, 419 (1981). 

[47] M. Christensson and M. Hindmarsh, Phys.Rev.D 60, 063001 (1999); M. Hindmarsh, 
these proceedings. 

[48] A. Berera and D. Hochberg, |cond-mat / 1 03447 . 

[49] D. N. Schramm, R. V. Wagoner, Ann. Rev. Nucl. Part. Sci. 27, 37 (1977). 
[50] S. Sarkar, Rept. Prog. Phys. 59, 1493 (1996). 

[51] V. A. Rubakov and M. E. Shaposhnikov Usp.Fiz.Nauk 166,493 (1996) [Phys.Usp. 
39, 461 (1996)]. 

[52] K. Jedamzik and G. M. Fuller, Astrophys. J. 423, 33 (1994); H. Kurki-Suonio K. 
Jedamzik and G.J. Mathews, Astrophys. J. 479, 31 (1997); K. Kainulainen, H. Kurki- 
Suonio, E. Sihvola, Phys.Rev.D 59 083505 (1999). 

[53] K. S. Thorne, Astrophys. J. 148, 51 (1967); S. W. Hawking and R. J. Tayler, Nature 
309, 1278 (1966); J. D. Barrow, Mon. Not. R. Astron. Soc. 175, 359 (1976). 

[54] G. Steigman, Ann. Rev. Astron. Astroph. 14, 339 (1976). 

[55] H. Kurki-Suonio and E. Sihvola, Phys. Rev. Lett. 84 3756 (2000); Phys.Rev.D 62, 
103508 (2000). 



31 



[56] J. B. Rehm and K. Jedamzik, Phys.Rev.Lett. 81, 3307 (1998); Phys.Rev.D 64, 
023510 (2001) 



[57; 

[58 
[59 
[60 
[61 
[62 

[63 

[64 
[65 

[66 
[67 
[68 
[69 
[70 
[71 
[72 
[73 
[74 
[75 



J. D. Barrow, Phys. Rev. D 35, 1805 (1987). 

E. W. Kolb, M. J. Perry, and T.P. Walker, Phys. Rev. D 33, 869 (1986); F. S. 
Accetta, L. M. Krauss, and P. Romanelli, Phys. Lett. B 248, 146 (1990). 

T. Appelquist, A. Chodos, and P. G. O. Freund, Modern Kaluza Klein Theories 
(Addison- Wesley, Redwood City, CA, 1987). 

M. B. Green, J. H. Schwartz, and E. Witten, Superstring Theory (Cambridge Uni- 
versity Press, Cambridge, England, 1987). 

J.K. Webb, M.T. Murphy, V.V. Flambaum, V.A. Dzuba, J.D. Barrow, C.W. 
Churchill, J.X. Prochaska, and A.M. Wolfe, Phys. Rev Lett. 87, 091301 (2001). 

M.T. Murphy, J.K. Webb, V.V. Flambaum, V.A. Dzuba, C.W. Churchill, J.X. 
Prochaska, J.D. Barrow, and A.M. Wolfe, Mon. Not. Roy. Astron. Soc. 327, 1208 
(2001). 

J. K. Webb, V. V. Flambaum, C. W. Churchill, M. J. Drinkwater, and J. D. Barrow, 
Phys. Rev. Lett. 82, 884 (1999). 

M. Giovannini, Phys.Rev.D 55 595 (1997). 

M. Gasperini, M. Giovannini and G. Veneziano, Phys. Rev. Lett. 75, 3796 (1995); 
Phys. Rev. D 52, 6651 (1995); M. Giovannini, Phys.Rev.D 59, 123518 (1999). 

M. Giovannini, Phys.Rev.D 56, 631 (1997). 

M. Giovannini, Phys.Rev.D 62 123505 (2000). 

M. Giovannini, Phys.Rev. D 64, 061301 (2001); 

M. Giovannini, |hep-ph/0104214l . 

P. A. M. Dirac, Nature 139, 323 (1937); Proc. R. Soc. London A 165, 199 (1938). 

P. Jordan, Z. Phys. 157, 112 (1959). 

E. Teller, Phys. Rev. 73 , 801 (1948). 

J. D. Bekenstein, Phys. Rev. D 25, 1527 (1982). 

Smith et al., Phys. Rev. D 61, 022001 (2000); J. C. Long et al., |hep-ph/0009062 
C. D. Hoyle et al., Phys. Rev. Lett. 86, 1418 (2001). 



32 



[76] G. German and G. G. Ross, Phys. Lett. B 172 305 (1986). 

[77] A. S. Goncharov, A. D. Linde and M. Vysotsky, Phys. Lett. B 147 279 (1993). 

[78] M. Gasperini, Phys. Lett. B 327, 214 (1994). 

[79] V. A. Kuzmin, V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B 155, 36 (1985); 
V. A. Robakov and M. E. Shaposhnikov, Phys. Usp. 39, 461 (1996) [Usp.Fiz.Nauk 
166, 493 (1996)]. 

^0] M. E. Shaposhnikov, JETP Lett. 44, 465 (1986). 

^1] M. Giovannini and M. Shaposhnikov, Phys. Rev. D 62, 103512 (2000). 

52] G. Greenstein, Nature 223, 938 (1969). 

^3] J. J. Matese and R. F. O' Connel, Astrophys. J. 160, 451 (1970). 

M] B. Cheng, D. N. Schramm and J. Truran, Phys. Rev. D 45, 5006 (1994); B. Cheng, 
A. Ohnto, D. N. Schramm and J. Truran, Phys. Rev. D 54, 4174 (1996). 

55] P. Kernan, G. Starkman and T. Vachaspati, Phys. Rev. D 54, 7207 (1996). 

56] D. Grasso and H. Rubinstein, Astropart. Phys. 3, 95 (1995) and Phys. Lett. B 379, 
73 (1996). 

57] K. Jedamzik, V. Katahnic, and A. V. Ohnto, Phys. Rev. Lett. 85 700, (2000). 

^8] J. Adams, U. H. Danielsson, D. Grasso, and Hector Rubinstein, Phys. Lett. B 388, 
253 (1996). 

[89] R. Durrer, M. Gasperini, M. Sakellariadou, and G. Veneziano, Phys. Rev.D 59, 
043511 (1999). 

[90] M. Giovannini, Phys. Rev. D 56, 3198 (1997). 

[91] M. Rees, Astrophys. J. 153, LI (1968). 

[92] P. M. Lubin and G. F. Smoot, Astrophys. J. 245, 1 (1981); P. M. Lubin, P. Melese 
and G. F. Smoot, Astrophys. J. Lett. 273, L51 (1983). 

[93] M. Bersanelh et. al. Esa D/Sci(96)3; home page of the Cambridge Planck Analysis 
Center (CPAC) ||http: / / www. mrao. cam. ac.uk/ projects / cpac /| . 

[94] E. Milaneschi and R. Fabbri, Astron. Astrophys. 151, 7 (1985). 

[95] D. Coulson, R. Crittenden and N. G. Turok, Phys. Rev. D 52, 5402 (1995); Phys. 
Rev. Lett. 73, 2390 (1994). 

33 



[96] A. Kosowsky and A. Loeb, Astrophys. J 461, 1 (1996). 

[97] D. Harari, J. Hayward and M. Zaldarriaga, Phys. Rev. D 55 (1997). 

[98] M.Giovannini, Phys.Rev.D 61, 063502 (2000). 

[99] M. Giovannini, Phys.Rev.D 61, 063004 (2000). 



[100 
[101 
[102 
[103 
[104 

[105 
[106 
[107 
[108 
[109 
[110 
[111 

[112 
[113 

[114; 

[115 



[116 
[117 



R. Jackiw and S.-Y. Pi, Phys.Rev. D 61, 105015 (2000). 

R. Jackiw, V.P. Nair, S-Y. Pi, Phys.Rev.D 62, 085018 (2000). 

C. Adam, B. Muratori, and C. Nash Phys.Rev.D 61, 105018 (2000). 

C. Adam, B. Muratori, and C. Nash, Phys.Rev.D 62, 105027 (2000). 

K. Kajantie,M. Laine, J. Peisa, K. Rummukainen, M. Shaposhnikov, Nucl.Phys.B 
544, 357 (1999). 

P. Elmfors, K. Enqvist, and K. Kainulainen, Phys.Lett.B 440 269, (1998). 

A. Vilenkin, Phys. Rev. D 22, 3067 (1980). 

V. Rubakov and A. Tavkhelidze, Phys. Lett. B 165, 109 (1985). 

V. Rubakov, Prog. Theor. Phys. 75, 366 (1986). 

A. N. Redhch and L. C. R. Wijewardhana, Phys. Rev. Lett. 54, 970 (1984). 

M. E. Shaposhnikov, Nucl. Phys. B 287, 757 {1987);ibid. 299, 797 (1988). 

K. Kajantie, M. Laine, K. Rummukainen and M. Shaposhnikov, Nucl. Phys. B 495, 
413 (1997). 

M. Carena, M. Quiros and C. Wagner, Phys. Lett. B 380, 81 (1996). 

M. Laine Nucl. Phys. B 481, 43 (1996). 

J. Cline and K. Kainulainen, Nucl. Phys. B 482, 73 (1996). 

B. Campbell, S. Davidson, J. Elhs and K. Olive, Phys. Lett. 297B , 118, (1992); 
L. E. Ibanez and F. Quevedo, Phys. Lett., B283, 261, 1992; J. M. Cline, K. Kain- 
ulainen and K. A. Ohve, Phys. Rev. Lett. 71, 2372 (1993); Phys. Rev. D 49, 6394 
(1993). 

M. Giovannini, Phys. Rev.D 59, 063503 (1999). 
M. Giovannini, Phys. Rev.D 58, 124027 (1998). 

34 



[118] D. Deryagin, D. Grigoriev, V. Rubakov and M. Sazhin, Mod. Phys. Lett. A 11, 593 
(1986). 

[119] M. Giovannini, Phys.Rev. D 60, 123511 (1999). 

[120] D. Babusci and M. Giovannini, Int.J.Mod.Phys. D 10, 477 (2001); Class. Quant. 
Grav. 17, 2621 (2000). 



35 



